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NAVAL ARCHITECTURE 


Location 
of ship rolling axis 


Ludwik Balcer 
Gdańsk University of Technology 


ABSTRACT 


In the paper is presented a definition of ship rolling axis and a method for determination of 
location of the axis. Location of rolling axis depends not only on ship mass distribution but 


E 


* | also on the so-called „added masses” which are equivalents of the hydrodynamic forces 
4 acting on ship hull during its motion. It has been assumed that the rolling axis is fixed with 
y 


respect to ship's hull, and its location is determined with accounting for the added masses 
accompanying ship rolling, whereas the influence of damping forces is neglected. For 
a representative group of ships appropriate coefficients of hydrodynamic forces were cal- 
culated on the basis of the experimental test results available in the subject-matter literatu- 
re. It has been assumed that the rolling axis passes through a common centre of ship's mass and relevant 
added masses. As a result proposed are simple approximate formulas which make it possible to calculate 
location of ship rolling axis on the basis of the typical ship data available at the preliminary design stage. 


Key words : ship seaworthiness, ship safety, ship hydromechanics, stability of floating units 


INTRODUCTION 


The ship rolling axis is a conventional notion which serves 
mainly to simplified representation of ship motions in the case 
when ship rolling prevails. Under the assumption that the loca- 
tion of ship rolling axis is known it is easy to determine trans- 
lations and accelerations in any ship's point, as then they are 
functions of distance of a given point from the axis in question. 
In some of the points the placing of some cargoes, measuring 
instruments, passenger cabins should be avoided or — in rever- 
se — they may appear most favourable for installing given de- 
vices. In order to estimate rolling axis location a ship physical 
model of two degrees of freedom (dof ) was assumed. How- 
ever if the location is assumed permanent the model can be 
taken as that of one dof only. In this model it is assumed that 
the considered ship rolls in calm water due to an external exci- 
tation free moment of appropriately selected parameters; it can 
also roll freely. In this paper the method of determining loca- 
tion of the rolling axis is described and — on this basis — sub- 
mitted are simplified formulas for calculation of its location, 
applicable to typical merchant ships. 


DEFINITION OF SHIP ROLLING AXIS 


It is usually assumed that the ship rolling axis coincides 
with X-axis of the coordinate system whose origin is fixed 
with respect to the ship in its mass centre [3,5]. It is also assu- 
med that the axis which lies on the ship plane of symmetry 
and initially is parallel to ship waterline, can be considered as 
the ship's main axis of inertia. With respect to it, calculated is 
the ship mass moment of inertia Ix, an element appearing in 


equations of motion. The axis is assumed motionless with re- 
spect to both the ship and space. The approximate character of 
the calculations makes it possible to use values of the inertia 
moments of the assumed form, i.e. (1, + Ixx), in spite of the 
fact, that the inertia moments are not corrected with taking 
into account a real location of rolling axis. It was this way 
assumed because the magnitudes of the added mass inertia 
moments of water, Ixx, are often determined on the basis of 
the measured free-rolling period of ship, hence the correct lo- 
cation of rolling axis has been already accounted for. This also 
concerns the approximate calculation formulas for both Ix and 
I; values. 


LOCATION OF SHIP ROLLING AXIS 


In order to determine location of ship rolling axis in accor- 
dance with the above specified assumptions it is postulated that 
the ship oscillates in calm water, having two degrees of free- 
dom : roll and sway. In this case the line on the ship's plane of 
symmetry, where lateral translations in the motionless, earth- 
-fixed coordinate system are equal to zero, can be considered 
as the rolling axis. The similar approach is presented in [2 , 7]. 


The following equation of ship's rolling 
is taken for considerations : 


т ўе+т»„ў+2н„ў +т„Ф+2н„Ф=0 
and after changing notation : (1) 


my tm, y 42u,, y +m, oP +2p,,0 =0 
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where : 


yg -lateral translation of ship mass centre 
y -lateral translation of the point O (Fig.1) 
Ф -ship rolling angle 
m  -ship mass 
Myy -added mass of water for ship sway 
туф - added mass of water for ship sway 
in result of rolling with respect to the point Og 
Шуу - damping coefficient for ship sway 
Hy - damping coefficient for ship sway 
due to rolling with respect to the point Oo 


This is the equation which represents the sum of projec- 
tions — on Y-axis — of the forces acting upon the ship, due to 
sway. In the equation (1) there are no expressions for restoring 
forces which really do not appear during sway. 


Fig. I. The coordinate systems and location of rolling axis 


In Fig.l are presented : the ship in a heeling position, the 
motionless, earth-fixed coordinate system Oo, Y ,Zo, and the 
movable, ship-fixed coordinate system O,X, Y. In the initial po- 
sition of the ship both the coordinate systems coincide. 

In Fig.1 the point O; represents the trace of the rolling axis 
and in accordance with the assumed definition it always rema- 
ins on Zo- axis of the motionless coordinate system. If in the 
ship-fixed reference system the quantity ZGw stands for the di- 
stance of the mass centre G from the point O, i.e. from the 
initial ship waterline, then the trace of the rolling axis, O,, splits 
the segment OG into two parts : 


ZGw= aw + bw 
For small heeling angles it can be written as follows : 


and у=аФ 


yG = - by® 


The quantities y determine shifts of the origin O of the 
movable coordinate system. The added mass туф and the dam- 
ping coefficients руу , руф are determined with respect to the 
origin О. The here considered damping coefficients цуу and 
Uyg take into account only dissipation of the energy generated 
by the radiation forces with neglecting the viscosity forces. 

As the segment ay = OO,, and by = GO,, hence ay< 0 
when the trace of the rolling axis, O,, lies over the initial wa- 
terline, i.e over the point O, and by < 0 when O, lies below the 
ship mass centre G (Fig.1). Now the equation (1) can be writ- 
ten as follows : 

-mb,o + my,,a,o + 
i » : (2) 
*2u,,a, P+M oP + 2u,4 =0 
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Knowing that ay = zGw - by , 
and introducing appropriate transformations one obtains : 


— by [(m+ my, )@ + 2u,,0 [= @) 


n (m, zo, emus )Ф -2 (uy zo Hye jò 

If the influence of damping is now neglected, i.e. руу = 0 
and Hyg = 0, then the following expressions for the distance 
bw of the rolling axis point O, from the ship mass centre G, and 
for the distance ay of the point O, from the initial waterline 
(i.e. from the point O, Fig.1) are obtained: 


MyyZGwt Myo 
by = 
M+Myy (4) 
ZGwIn — Myo 
ay = —— 


m+ Myy 


Therefore, with the lack of damping the rolling axis loca- 
tion is fixed for a given rolling period, and it accounts for the 
influence of the inertia forces resulting from the water surro- 
unding the ship, here represented by the added masses myy and 
туф. If the ship having the same rigidity, rolled in the „free 
space”, i.e. not being immersed in the water, then туу and туф 
would equal zero. In such case the rolling axis would go thro- 
ugh the mass centre С, and then it would be : 


by = 0 and ay = ZGw , respectively. 


In order to determine rolling axis location for a considered 
ship it is necessary to calculate the added masses myy and туф. 
For the quantities it is hard to find approximate value, which is 
not the case for e.g. rolling or heaving. The quantities can be 
calculated analytically by using the ,,strip theory” and apply- 
ing the Lewis frame-like forms [3, 5, 7]. However this requires 
using a special software together with many initial data of 
a considered ship. 

For these reasons a simplified method is proposed below. 
It can be useful for the approximate determination of rolling 
axis location on the basis of the ship's data usually available 
for a designed ship. 


A METHOD OF DETERMINING 
LOCATION OF ROLLING AXIS 


In order to determine the location of rolling axis in which 
an influence of added water mass is accounted for, a simplified 
ship representing a given one has been assumed. The simpli- 
fied ship is a cylinder of rectangular frame-section and its pa- 
rameters are so selected as to obtain its roll and sway motions 
as close as possible to those of the basic ship. Hence the ship is 
of the same mass m, B/d ratio and free-rolling period t. The 
centre of gravity G of the simplified ship is located, proportio- 
nally to its main dimensions, in the same point as that of the 
real ship. It has also the same transverse moment of inertia of 
waterplane area, which indirectly gives the same small meta- 
centric radiuses. The last condition makes it possible to deter- 
mine the main dimensions, Lg, Bg, ds, of the simplified ship on 
the basis of the real ship's parameters. It was further assumed 
that a change of location of rolling axis of the simplified ship, 
resulting from added mass influence, correctly describes the 
corresponding change in the basic ship. 

Values of the added masses and damping coefficients for 
the ship of rectangular frame sections (as well as for some other 
simple forms) are published in [8] and [3] in the form diagrams 
of the dimensionless values m'yy, Шуф, of the added masses 
Myy and myg, in function of the dimensionless ship rolling fre- 
quencies 0)". 


In order to determine m'yy it is also necessary 
to know B/d ratio. 


Now it is possible to describe successive steps for determi- 
nation of ship's rolling axis location in compliance with the 
above made assumption. 


Step 1. 
Selection of ship's data 


The demanded ship's data are as follows : 
Ly, -length between perpendiculars [m] 
B -breadth [m] 

d -draught [m] 

H -depth to upper deck [m] 

m -ship mass [t] 

d -block coefficient of ship's hull [-] 

a - waterplane coefficient [-] 

ZG -height of the ship mass centre 
over the ship base plane [m] 

GM - initial transverse metacentric height 
of the ship, not corrected 
for free- surface influence [m] 


In order to obtain calculation results of a sufficiently gene- 
ral character the data of different ships in various loading con- 
ditions were used (Tab.1). 


Step 2. 
Determination (calculation) of ship's rolling period c 


The best solution is to have a value of the rolling period 
measured on the real ship. However such data are hard availa- 
ble. An approximate value of ship's free-rolling period can be 
calculated by using the following relationship : 


1„ +1 
= e (5) 
mgGM 
To determine values of I, and I, x 
the expressions given in [6] were used : 
m 
aC (в2+422) L5 0.31, (6) 


In further considerations the quantity (I, + Ixx) is deemed 
a sufficient approximation of the total moment of inertia of 
ship mass and added mass of water, respective to ship rolling 
axis, as it was postulated in its definition. 


Step 3. 
Calculation of the main dimensions L, , B, , d; 
of the simplified ship 


The ship breadth B, is determined by assuming that values 
of the moment of inertia, Iwx, of waterplane areas of the real 
ship and the simplified one are equal. The Iyx of the real ship 
waterplane is determined by approximating ordinates of its 
contour, yw(x), by means of the parabola whose order is expres- 
sed by the waterplane coefficient о, [4]: 


a 


B X -« 
y= Taa 


2 (+о)(ї+20) 


Basing on Iwx given by (7) 
and in accordance with the following relationships : 


В В 


Ela Ls Bsd =LBdôð =. 
12 d, d 


one obtains the expression for the breadth Bs, draught ds 
and length L, of the simplified ship : 


6B o 


8 (i+a)(1+2a) 


L, =| —-| L8 
B B, 


s 


Step 4 and 5. 
Calculation of œ', read-out of m'y, and m'yo , 
and calculation of my, and myo 


The form of the quantities to be calculated, in compliance 
with the previously applied symbols, is as follows : 


B m 
@'= 0) 3 m'y = УУ 
2g pA, 
А П уф; 
= ———— (9) 
Y" pA,B, 


my,-L, my; Myo = Ls + Myq; 


The values туу; [t/m] and туф; [tm/m] concern the ship of 
the one-meter length and of the assumed, simplified frame 
section. 


The remaining symbols stand for : 


p -— water mass density [t/m?] 

g -— acceleration of gravity [m/s] 

B, — breadth of the simplified ship [m] 

Ag — ship's frame area [m] (here : A; = Bg dg) 


Having read values m'yy and m'yF according to w', one calcu- 
lates, values туу; and туру, and next туу and myr, by using (9). 


Step 6. 
Calculation of rolling axis location 


The calculations were carried out for the simplified ships. 
First, by using (4) the segments ays and bys were calculated 
and then applied in their relative form by making use of the 
following relations: 


to obtain: aw = (aws/Bs)‘B and by = (bws/Bs)B 


Next, the values of the added masses, Myy and туф, Seg- 
ments aws and bws, as well as ratios ay,/B, and bws/B; were 
calculated on the basis of the above presented assumptions and 
formulas. 

The calculations were performed for 9 ships in different 
loading conditions. Tab.1 contains the input data of the ships, 
and the calculation results of the relative values of rolling axis 
location are shown in Tab.2 ( according to the increasing sequ- 
ence of ZGy/B values). 

The segment ay = OO, was selected as the quantity determi- 
ning location of rolling axis. The segment ay, being the distance 
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Tab. 1. Input data of 9 selected ships in different loading conditions 


Saeed Maan ie 
0.672 
iso | 1940 | 36] ssa2 [oms оз [nase] ia 
TE 

ful 
63 
66 


ac 134.00 | 21.40 19971 | 0.800 


Jerryship 
STENA 
Jerryship 


ROGALIN empty 


empty 


ship 
WARSZAWA 
dry cargo ship 
MAESTRO 


tanker 


GIEWONT II nd 


ballas 


bulk carrier 
UNIWERSYTET 


grain 


ballast 182.00 
Ro-Ro ship 


POZNAN 182.00 
"ne 
ishing trawl 


of rolling axis from waterline, correctly describes the influence 
of the added masses determined relative to the point O, (Fig.1), 
on location of rolling axis. 

It was also assumed that the relative quantity aw/B was most 
tightly connected with the relative distance of the ship mass 
centre from the same waterline, zGy/B. 

The relationship ay/B = f(zGw/B) is shown in Fig.2. This 
diagram reveals good correlation between the selected quan- 
tities. A significant part of the irregularities observed in the 
diagram can be deemed as a result of the accepted simplifying 
assumptions, as well as inaccuracy in reading-out intermedia- 


“0.4 


-0.3 -0.2 -0.1 0.0 0.1 0.2 


Fig. 2. Location of ship rolling axis respective to ship waterline 
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ful | 272.00 | 0.826 | 455 


272.00 | 43.40 | 10.45 | 102540 | 0.851 | 0.807 | - 0.28 
full 205.00 | 30.50 | 12.09 | 62450 | 0.871 
205.00 | 30.50 35620 | 0.818 


Tab. 2. Calculation results of rolling axis 
location for the selected ships 


D 


0.087 - 0.087 


038 
043 


0.140 | -0.062 


te values. Though the diagram has been obtained for the sim- 
plified ships, it can be considered as that correctly describing 
the selected relationships also for the real ships. The straight 
line described as : 


а,/В = 0.432 (zGw /B) + 0.102 
fits into the relationship aw/B = f(zGw/B) the best. 


ZGw/B 


| - 6.92 | | 15.6 | 
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©! © 
© 
чә 
чә 


© 
оо 
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N 


© 
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.087 


0.105 


(10) 


After a simple analysis of (10) it can be demonstrated that 
the rolling axis is located approximately 0.1 B under the water 
when the ship mass centre lies on its surface, and that the rol- 
ling axis is located on the water surface when the ship mass 
centre is abt. 0.24B over the water, and finally, that the rolling 
axis goes through the ship mass centre when it is located abt. 
0.18B under the water surface. To make the diagram more cle- 
ar, also the straight line ZGw/B = f(zgy/B) was drawn in it. 

By taking into account that zGy = d - zg the following ap- 
proximate expressions for the distance of ship rolling axis from 
ship's waterline and ship mass centre, respectively, can be ob- 
tained : 


ay = 0.43 (d - zg) + 0.1B 
(11) 


or 


by = 0.57 (d - za) - 0.1B 


This is a simplified form of the expressions (4) in which 
the hydrodynamic coefficients my, and туф functionally de- 
pend also on the rolling frequency o. The, observed in Fig.2, 
small scatter of the points, 1.e.values determined by (4), indi- 
cates that rolling axis location weakly depends on other para- 


meters than those represented by the coordinates shown in the 
diagram in question. In the obtained formulas the dependence 
on the rolling frequency has place indirectly only because the 
calculations for all the selected ships have been carried out with 
taking into account their free-rolling frequencies, hence the 
diagram of Fig.2 and the formulas (10) and (11) obtained on its 
basis, concerns just such frequencies. 


FINAL REMARKS AND CONCLUSIONS 


O Ship rolling motions are of a great importance in conside- 
ring the problems of ship stability, strength, cargo lashing, 
comfort and safety of passengers and crew members; they 
are also important for correct operation of ship equipment. 

O For considering and solving many of the above mentioned 
problems it is usually sufficient to assume the ship to be an 
object of one degree of freedom. When such approach is to 
be used correct determination of location of ship rolling 
axis is always crucial. 

O Due to the water supporting the ship during its motions, 
which forms constraints of a special type, it is possible to 
assume the free ship rotation axis to be a straight line pas- 
sing through — or close to — the ship mass centre. However 
the total ship mass consists of the mass of the ship itself 
and the supplementary, „added masses” of water, associa- 
ted with relevant motions. These conventional masses af- 
fect location of ship rolling axis in accordance with places 
of their occurence. 

O Ifthe ship rolling axis is so defined as it has been propo- 
sed in this paper the usually assumed location of rolling 
axis can be improved; the corrections might reach even 
0.17 0.2B as it was shown in this paper. 

O It is also important that the proposed formulas, (10) and 
(11), concern resonant rolling, as in real conditions the ship 
responds to instantaneous excitations with motions of na- 
tural frequency, and in irregular waves it ,,selects" — from 
an encountered wave spectrum — mainly the frequencies 
close to natural frequencies of its own motions. 

O The proposed simple approximate formulas make it possi- 
ble to calculate location of ship's rolling axis only on the 
basis of the typical ship data available at the preliminary 
design stage. 


NOMENCLATURE 

aw - distance from rolling axis to waterline of a real ship [m] 

aws - distance from rolling axis to waterline of a simplified 
ship [m] 

by - distance from rolling axis to mass centre of a real ship [m] 

bws - distance from rolling axis to mass centre of a simplified 
ship [m] 

dof - degress of freedom 

D - ship displacement [T or kN] 

g - acceleration of gravity [m/s?] 

GM - initial metacentric height ofa ship, not corrected for free- 
surface influence [m] 

Ix - ship mass moment of inertia relative to X- axis [tm?] 

Lx - added mass moment of inertia relative to X- axis [tm^] 

L, B, d, H - main dimensions of a ship [m] 

m - ship mass [t] 


m'yy, m'yo - dimensionless coefficients of added mass of water [-] 

Myyj, туфј- unitary coefficients of added mass of water [t/m, tm/m] 

Myy, туф - coefficients of added mass of water [t, tm] 

5 - index of a simplified ship 

YG Y - swaying translations of ship mass centre G and point 0, 
respectively [m] 

ZG - height of ship mass centre over base plane [m] 

ZGW - height of ship mass centre over waterplane [m] 


[»/ - waterplane coefficient [-] 

ô - hull block coefficient [-] 

T - period of ship's free rolling [s] 

o - heeling angle [rad or deg] 

[0] - angular frequency of ship's free rolling [1/s] 
[ny - dimensionless ship rolling frequencies [-] 
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On 11 December 2003 at Faculty of Ocean Engineering 
and Ship Technology, Gdańsk University of Technology, 
held was the last-in-that-year seminar of the Regional 
Group of the Section on Exploitation Foundations, Ma- 
chine Building Committee, Polish Academy of Sciences. 


During the seminar three papers were presented : 


e A global model of environmental safety of ship power 
plant — by R. Liberacki (Gdansk University of Tech- 
nology) 

e Research on water-lubricated bearings of ship main 
shafting — by W. Litwin (Gdansk University of Tech- 
nology) 

e Anattempt to diagnosing the main propulsion systems 
of the ferry ships "Polonia" and "Pomerania" by me- 
ans of the FAM-C method — by A. Gebura (Air Force 
Institute of Technology, Warsaw) 


The seminar was a good occasion to offer wishes to 
Prof. J. Lewitowicz in his 70” birthday and to pass gratu- 
lations to him due to his outstanding achievements in the 
area of exploitation of technical objects, especially in aero- 
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On a risk-based method for safety 
assessment of a ship in critical conditions 
at the preliminary design stage 


Current report 


Mirostaw Gerigk 
Gdansk University of Technology 


ABSTRACT 


The paper presents some current results of the research on development of a method for 
safety assessment of ships in critical conditions at the preliminary design stage. The me- 
thod is expected to make it possible to perform the multi-objective and multi-parametric 
investigations to achieve the required and optimum levels of safety by using the risk-based 
and formal approaches to ship safety. During the design analysis influence of parameters 
associated with the hull form, arrangement of internal spaces, loading conditions, position 
and extent of damage, cargo shift and weather impacts on safety of ships in critical condi- 


tions will be taken into account. Current capabilities of the method are illustrated by several calculation 
examples of influence of different arrangements of internal spaces of a cargo ship on values of the IMO 
formal subdivision indices such as the "Attained Subdivision Index A” and "Local Subdivision Indices AA; ". 
The method is still under development and its final structure will be presented in the near future. 


Key words : ship safety, formal safety assessment, design for safety, ship hydromechanics, survivability 


INTRODUCTION 


During the last decade many European institutions conduc- 
ted the investigations related to the ship safety problems [1, 2, 3, 
4, 5]. Recently, the HARDER research programme has been set 
up to revise the SOLAS Convention, (Chapter II-1 Parts A, B 
and B-1), it has been aimed at solving the following problems [8] : 


» reviewing model tests and methods concerning safety eva- 
luation and assurance including the development of the 
methods contained in A.265 IMO document [14], and also 
that applicable to Ro-Ro ships 

> survival prediction for various ships in waves including 
observed mechanisms of capsize 

» estimation of probability of sea state occurrence at the time 
of a casualty 

> development of the so called survival factor “5”. 


Also, at the Chair of Ship Hydromechanics, Faculty of 
Ocean Engineering and Ship Technology, Gdansk University 
of Technology two projects have been conducted regarding 
a new method and model for the safety estimation of ship in 
critical conditions, as well as a model for direct risk assessment 
of ships in critical conditions [6 , 7]. The major research issues 
associated with the above mentioned projects were : 


** damage stability modelling methods 

** large scale flooding 

* dynamic effects due to internal (ballast and/or cargo shift) 
and external (waves, wind) impacts 

* development of survival criteria for the ships in damage 

condition 
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* investigations on the safety assessment 
* and example preliminary designs. 


The critical conditions have been defined as the ingress of 
external water into any group of watertight compartments of 
a ship, associated with both internal (cargo and/or ballast shift) 
and external (waves and wind action) impacts. 

The knowledge base for the research includes: naval archi- 
tecture, ship hydromechanics, system approach to safety, safe- 
ty case/formal safety assessment (FSA) methods, IMO regula- 
tions for ship safety. 

This paper is devoted to presentation of some current re- 
sults of the conducted research projects regarding development 
of a method for safety assessment of ships in critical condi- 
tions at the preliminary design stage. 


MAJOR ASSUMPTIONS FOR THE 
MODERN APPROACH TO SHIP SAFETY 
IN CRITICAL CONDITIONS 


The modern approach to ship safety in critical conditions 
is based on the general assumption as follows : 


The first assumption 
The factors affecting the ship safety in critical conditions 
may result from different sources : design, operation and ma- 
nagement. 


The second assumption 
The factors affecting the ship safety may come 
from different groups of factors such as : 
* Ist-the external factors not related to the environment : 


general human factor, control systems, technical me- 
ans, legislative actions. 
Ф 2nd - the environmental factors including the wind, wa- 
ves and current. 
* 3rd - connected with the ship : 
= ship including hull, propeller and rudder 
= cargo including arrangement of internal spaces, cargo 
and ballast distribution and loading condition 
= inter-related parameters and characteristics directly 
connnected with ship safety: intact stability, dama- 
ge stability, dynamical stability in damage condi- 
tion, and survivability. 
* 4th - connected with ship operation: the ship manage- 
ment system including the ship safety system. 
Ф Sth-connected with the shipping company’s safety ma- 
nagement system and safety culture. 
* 6th -the human factor including both the psychological 
and physical predispositions, character, morale, integrity, 
knowledge, experience and training degree. 


In the current work the factors ofthe first four groups have 
been taken into account in the method for assessment of ship 
safety in critical conditions. 


The third assumption 
The modern approach to ship safety is connected with com- 
bining: the system approach and Formal Safety Assessment 
(FSA) method which contains : hazard identification, risk as- 
sessment and risk reduction. 


By taking into account these elements the method in ques- 
tion has gradually been developed to include the following 
components (shown also in Fig.1) : 


design requirements, criteria and constraints 
risk acceptance criteria 
safety objectives 
ship and environment definition 
design analysis : 
= hazard identification 
= hazard assessment 
= scenario development 
= hydromechanical design analysis 
(intact stability, damage stability, 
dynamical stability in damaged condition) 
risk assessment 
risk acceptance criteria 
risk reduction (mitigation measures) 
modification of design 
ship safety assessment 
decisions on ships safety. 


For estimating the risks the following techniques 
have been selected to be used within the method : 


= ALARP (As Low As Reasonable Possible) concept 
c F-N curve: the frequency of fatalities 

in function of number of fatalities 
= Risk acceptance matrix. 


For the hazard identification the hazard and operability (HA- 
ZOP) studies have been used as the most effective technique. 
The event tree analysis (ETA) has been used for the scenario 
identification. 

The statistics and investigations into serious casualties, 
documented in [9], have been assumed the major source of 
information on hazards and risks involved in shipping. It has 
been further assumed that the design method should combine 
both global approach and technical approach [13]. 


containing the elements just above enumerated. The technical 
approach has been connected with the developing of : 


The global approach has concerned the method framework 


logical structure of design system 

design requirements, criteria and constraints 
logical structure of computational model 
analytical and numerical methods 

application methods regarding the intact stability, 
damage stability, dynamical stability in critical 
and damage condition. 
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For risk management considerations the top-down risk 


management method has been applied which is suitable for 
design for safety at the preliminary design stage. 


METHOD OF SHIP DESIGN FOR SAFETY 
IN CRITICAL CONDITIONS 


A worked-out safety estimation method for ships in critical 


conditions is associated with solving a few problems regar- 
ding the naval architecture, ship hydromechanics and ships 
safety, and it is novel to some extent. 


The logical structure of the method and design system com- 


bining the global approach and technical approach is presen- 
ted in Fig.1. 


Risk Acceptance Criteria 


Safety Objectives 


Cost/Benefit Constraints 


Main Design Requirements 


Design Criteria & Constraints 


Ship & Environment Definition 
Hazard Identification 
& Hazard Assessment 
Design Analysis Scenario Identification 


Risk Reduction Hydromechanic 
Design Modification 


Design Analysis 
N 


[2] 


Risk Control 
Risk Tolerable 


Risk Assessment 


Decisions on Safety 


Fig.1. Logical structure of the design system and computational 
model combining the global approach and technical approach 


According to Fig.1, the ship and environment are defined 


as objects of hydromechanics, described by a set of parame- 
ters. The "Hydromechanic Design Analysis" module contains 
the design methods in which both the functions and procedu- 
res associated with solving particular ship hydromechanics pro- 
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blems are used. The "Risk Assessment" module includes the 
methods which combine both the "hydromechanics" and "risk 
assessment" functions and procedures. The "Main Design Re- 
quirements" module is intended to be consisted of : 


4 


* general requirements 

IMO regulations 

requirements of classification societies and 
requirements of conventions 

(SOLAS'90 and SOLAS' 95. 


X 


= 


X 


= 


X 


S 


The current set of requirements included into the database 
has so far concerned the IMO regulations only. The database 
will eventually include both the risk acceptance criteria, safe- 
ty objectives, main requirements and design criteria and con- 
straints. 

By using the method 
two separate design processes can be initiated: 


$ iterative approach 
$ parametric investigations. 


From the design process point 
of view the method consists of two sub-methods [10]: 


* parametric method - when intact stability and damage sta- 
bility characteristics are calculated 

* semi-probabilistic or probabilistic method - for solving the 
problems related to the survivability and risk assessment. 


The following design options of the computational model, 
based on the IMO regulations [15], are currently offered : 


ж calculation of ће "Attained Subdivision Index A" 
* calculation of the "Local Subdivision Indices AA;" 
** calculation of the "Probability of Oil Out-flow" 


Moreover it is possible to use the option : calculation ofthe 
"Probability of Capsizing in Critical Conditions", not covered 
by the IMO regulations, 1.e. the direct risk assessment method 
which is still under development. 


At the current stage of development of the method 
it has been assumed that the method : 


+ has scientific and research features 
rather than purely design one 

* isuseful for design of the conventional 
cargo vessels only. 


The method is still under development 
and its final structure will be presented in the near future. 


COMPUTATIONAL MODEL 
OF SHIP DESIGN FOR SAFETY 
IN CRITICAL CONDITIONS 


The most important features 
of the computational model in question are as follows : 


A itis open 

A it has a hybrid-modular structure 

A it has a common library of analytical 
and numerical methods 

А it has a common library 
of application methods. 


From the practical point of view the computational model 
is based on a dynamical database concept, and it is original to 
some extent. The database is modular and related to the logical 
structure of the computational model. The database enables to 
provide the safety estimation when ship’s hydromechanical cha- 
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racteristics are estimated by using : the numerical calculations 
(direct methods), model tests results, results from the full scale 
trials, combined empirical and numerical calculations (semi- 
-direct methods) or empirical calculations (indirect methods). 

The computational model has been prepared for the PC 
Windows-based workstations. The application software is pro- 
grammed in Fortran, Pascal and C++ languages. The detail 
structure of the data base including the input data, outputs, gra- 
phics and ship design data will be presented in the near future. 

At the current stage of work it is possible to use the me- 
thod and computational model for the research purposes only. 
A version for implementation for the research and design pur- 
poses would be available soon. 


SOME EXAMPLES 
OF CURRENT CAPABILITIES 
OF THE METHOD IN QUESTION 


Below presented are a few examples of using the method 
for calculation the probability of flooding any group of com- 
partments, optimisation of the "Attained Subdivision Index A" 
and calculation of the "Local Subdivision Indices AA;" in com- 
pliance with IMO relevant regulations. 


Example 1 


Calculation of the probability 
of flooding any group of compartments 


At the design stage the estimation of the probability of sur- 
vival of flooding any group of compartments is connected with 
calculation of the "Attained Subdivision Index A" [15]. The 
risk assessment is associated with satisfying the criterion : 


AR 
where : 
А = Pi Sj 


Both the indices, A and R, are calculated according to the 
formula accepted by IMO. In the presented example the for- 
mula given by the Resolution MSC 19/58 — "Subdivision and 
damage stability of cargo ships over 100 m" has been applied. 


A typical set of ship initial data 
for the given ship may look as follows : 


main parameters 

e stability data, hull form 

e arrangement of internal spaces 
(like the presented in Fig.2) 

e loading condition 

e set of permeabilities. 


Arrangement of internal spaces : Main Parameters : 


- Transverse Subdivision Considered Lg = 174.95 
Lbp = 163.00 
- Longitudinal Subdivision Considered B = 26.50 
- Horizontal Subdivision Considered H = 14.20 
d = 9.00 


ү. Gees 
aa 


- 86.29 


i ns === БЫ [ШР 
s qs fs fe rs |» 2 


- 86.29 


Fig.2. An example arrangement of internal spaces of a containership 


According to the performed calculations value 
of the "Attained Subdivision Index A" is as follows [16] : 


A= >ААр{ = 0.5501 
The Required Subdivision Index is: 


R = (0.002 + 0.0009 - Ls)! = 


= (0.002 + 0.0009 - 174.98)! = 0.54230 


where : 
Ls = 174.98 [m] 
The final result satisfies the relation A > R, as : 


0.5501 > 0.5423. 


Example 2 
Optimisation of the "Attained Subdivision Index A" 


For more advanced procedures of design for safety of 
a ship in damage condition the following approach can be ap- 
plied. In Fig.3 the arrangement of internal spaces of a cargo 
ship is presented. During the computer simulation three design 
versions regarding the number and position of transverse bulk- 
heads were taken into account. Tab.1 presents the results of 
optimisation of the "Attained Subdivision Index A" for all the 
design versions (Fig.3). During the optimisation process the 
preliminary positions of some bulkheads were significantly 
changed. 


Design version 1 


-78.60 77.10 


Fig. 3. Three different preliminary arrangements 
of internal spaces of the considered cargo ship 


The optimisation method consited in maximization of the 
objective function represented by the "Attained Subdivision 
Index A". Its aim concerned reaching possible maximum valu- 
es of the p; : s; factors for each group of the watertight com- 
partments. One of the constraints was that the s; value should 
never be equal to null. The major starting design conditions 
were : number and positions of bulkheads. 

The investigations showed that a larger number of bulkhe- 
ads not necessarily guarantee much higher values of the index 
"A". Moreover, by applying one bulkhead more the p; values 
become greater for the groups including two or more single 
compartments, but in the same time the conditional probabili- 
ties s; decrease. The values of p; and s; factors are interrelated 


Tab. 1. Results of optimisation of the "Attained Subdivision Index A" 


Ship type: general cargo 
Design criterion: optimisation 
of the "Attained Subdivision Index A" 


No. of iteration Index A value [-] 
I 0.6050 


but in an irregular manner. In general, one bulkhead more or 
less may only slightly affect the "Attained Subdivision Index 
A" values. The small changes in the positions of bulkheads do 
not give much different values of the index "A", but the local 
subdivision indices may be changed significantly. 


Design version 1 


Design version 2 


Example 3 
Calculation of the ,,Local Subdivision Indices AA "nu 


The calculation of the local subdivision indices is a kind of 
optimisation procedure which may provide the same level of 
safety for each of the considered watertight compartment. The 
local subdivision indices were calculated according to the fol- 
lowing formula [11]: 


AA; = (È pi si) / È pi) 


The optimisation process of the indices AA; is connected 
with the iterative shifting of the bulkheads against each other 
until the AA; values become more or less equal, possibly for all 
the compartments. In Tab.2 the local subdivision indices AA; 
are presented for four performed design iterations. The calcu- 
lations concerned the design version "1" of the general cargo 
ship presented in Fig. 3. 


-60.00 
-40.00 
0.4953 


-16.69 
15.00 
0.7457 


Tab.2. Results of the optimisation of the local subdivision indices AA; for the 
design version "A" is presented in Fig. 3 
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In Fig. 4 the graphical results 
of the computer simulation in question are presented. 


Local subdivision indices 


E 151 iteration 


0.914 


Relativ local 
subdivision 


1 2 3 4 5 6 7 
Compartments 


Local subdivision indices 


ш 2nd iteration 


Relativ local 
subdivision 


Compartments 


Local subdivision indices 


E 3rd iteration 


Relativ local 
subdivision 


Compartments 


Fig. 4. Results of the optimisation of the local subdivision indices AA; for 
the design version "1" of the general cargo ship, presented in Fig. 3 


RESEARCH PROJECTS UNDER WAY 


One of the projects concerns estimation of the probability 
of capsizing a ship in critical conditions by using the direct risk 
assessment method applied to stability problems. 

Another deals with consideration of a few novel safety as- 
sessment approaches with a view of implementing them into 
a safety system for ships in critical conditions to support the 
application of formal safety assessment. The following novel 
safety assessment and decision-aiding tools applicable to safe- 
ty management system are considered [14]: 


* expert judgements based on the Taguchi method 

** asafety based decision support system using either know- 
ledge-based expert systems or artificial neural network tech- 
niques 

* a fuzzy-logic-based synthesis incorporating the Dempster 
Shafer approach for making multiple-attribute decision 

* an integration of approximate reasoning approach and evi- 
dential reasoning method in design evaluation. 


SUMMARY 
Some results of the research projects on an integrated safe- 


ty estimation method and model for assessing the safety of 
a ship in critical conditions has been reported. The method is 
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directed towards the ship safety estimation in critical condi- 
tions at the preliminary stage of design. By using the method/ 
/model the hazard and risk assessment may be performed ac- 
cording to the IMO regulations for cargo ships. 

So far, the method and model have been used for investiga- 
ting the ship safety from the damage stability and survivability 
point of view. The method/model can be used for guiding ship 
subdivision for safety. Example results of the investigations of 
different arrangements of internal spaces of a cargo ship inclu- 
ding either transverse or combined subdivision, have been pre- 
sented. The damage stability, survivability and risk assessment 
calculations were carried out for each case. 


The method is still under development 
and its final structure will be published soon. 
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NOMENCLATURE 

A - attained subdivision index 

Ls - subdivision length 

Pi - probability of flooding the group of compartments 
under consideration 

R - required subdivision index 

Si - probability of survival after flooding the group of 


compartments under consideration 


AA; - local subdivision index related to each group of 
compartments under consideration 

АА, - components of the "Attained Subdivision Index A" for 
each group i of compartments 

X1, Хә - positions of bulkheads aft and fore, respectively 
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24" Symposium 
of Ship Power Plants 


On 27 and 28 November 2003 in Swinoujécie, a Po- 
lish port town on the Baltic Sea coast , was held 24th in- 
ternational meeting of the scientific workers engaged in 
the designing, testing and operating of ship power plant 
devices. The meeting was organized by the Institute of 
Technical Operation of Ship Power Plants, Maritime Uni- 
versity of Szczecin. 


The Symposium's agenda contained presentation 
of 44 papers split into 5 topical groups : 


= Design and operation of ship power plants (18 papers) 

= Renovation of ship power plants (2 papers) 

= Failures, tribology and diagnostics (10 papers) 

= Control of working processes and computer 
simulation (9 papers) 

= Solid/liquid separation (5 papers). 


The authors of the papers represented 7 Polish coastal 
scientific research centres and 3 Russian ones. Maritime 
University of Szczecin, Gdansk University of Technolo- 
gy and Technical University of Szczecin most contributed 
to preparation ofthe conference materials. The participants 
had the opportunity of being acquainted with current pro- 
duction achievements of such companies as : MAN-B&W, 
Caterpillar, Alfa Laval, Bosch Rexroth, Volvo Penta and 
Olympus, and of visiting floating objects of Polish Navy 
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16. Gerigk M. : A risk-based method for ship safety assessment at 


the preliminary design stage. Proceedings of 8th International 
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Miscellanea 


Agreement on cooperation 


In July 2003 the agreement on scientific technical coope- 
ration in the area of ship electric engineering between Polish 
Register of Shipping (PRS) and the Department of Ship Electric 
Engineering, Gdynia Maritime University, was signed by 
Mr. J. Jankowski, D.Sc., General Director of PRS, and Prof. J. 
Mindykowski, Head of the Department. 

Realization of the cooperation is intended on one hand to 
support developments in the field of ship electric power engi- 
neering, and on the other hand to help in developing and moder- 
nizing rules for construction and classification of ships and float- 
ing objects for sea and inland navigation. 


The common activity 
to be undertaken will concern : 


production and distribution of electric energy 

quality parameters of the energy 

influence of electric energy consumers on its quality 
influence of quality of electric energy on operation of 
its consumers 

instrumentation for measurement and control of elec- 
tric energy quality. 


Gug 


% 


Within the scope the cooperation 
will be realized by means of : 


reviewing topics proposed by each of the parties and 

common choice of priorities 

common carrying out measurements of electric ener- 

gy parameters on selected objects 

determining research directions on the basis of earlier 

obtained results 

common analyzing research results from the point of 

view of their possible application into PRS rules and 

instructions 

> analyzing possible improvements of ship electric plant 
control & measurement systems, resulting from rese- 
arch work 

> elaboration of a portable electric power analyzer ap- 

plicable to the new and existing ships not equipped 

with a modified permanent measurement & control sys- 

tem of electric power plant. 


v у у Y 


Finally common publishing activity in the range 
of the cooperation topics has been agreed. 


POLISH MARITIME RESEARCH, No 1/2004 13 


Miscellanea 


Marine propulsors - history and development 


Przemystaw Urbanski 


PEDNIKI 
OKRETOWE 


HISTORIA I ROZWOJ 


In the year 2001 the Shipbuilding & Shipping publishing 
house in Gdansk, edited (in Polish) the book on ,,Marine 
Propulsors - their history and evolution” written by Prze- 
myslaw Urbanski, retired Assist. Prof., D.Sc., M.E. The book 
resulted from the research undertaken by the author in co- 
operation with Hydrodynamics Department, Faculty of Oce- 
an Engineering and Ship Technology, Gdansk University of 
Technology, and Polish Nautological Association in the years 
1998 and 1999. In the this way elaborated monography the 
author has presented the entire image of history and deve- 
lopment of marine propulsors, understanding under this 
notion any device applicable to direct propelling various flo- 
ating objects beginning from the most primitive ones like 
rafts and boats to the most specialized e.g. hovercrafts or 
racing motor boats. Due to its popular-science character the 
book is dedicated to a wide spectrum of readers professio- 
nally or emotionally engaged in shipbuilding, yachting and 
water sports. 


Ressel's screw propeller (1812) 


The chronologically presented review of various techni- 
cal concepts and events enables to properly perceive genesis 
and ways of development ofthe contemporary propulsors as 
well as the role played by the gradually collected experience 
as well as by the development of engineering sciences. 


The entire scope of the book is contained 
in the eight chapters : 


I.  Punt-poles, Oars and Paddles 
П. Sails and Other Wind Propulsion Devices 
III. Paddle Wheels 
IV. Screw Propellers : 
e The origin, Design Types and Development of Screw 
Propellers Before 1836 
e F.P. Smith's and J. Ericsson’s Propellers. The Begining 
of the Era of the Screw-propulsion 
e The Screw-propulsion on the Sailing Ships 
e The Evolution of Screw Propellers after 1836 
V. Hydraulic Jet Propulsion 
VI. Vertical Axis Propellers — Cycloidal Propellers 
VII. The Idea of Nature — Inspired Propulsion Devices 
VIII. Past and Present Unconventional Propulsion Devices 


In the largest, 4th Chapter devoted to screw propellers, 
the geanalogical tree originally elaborated by the author, 
which illustrates paths of development of propulsors of the 
kind, is attention deserving. Roots of the tree go back to in- 
ventions of Archimedes, Leonardo da Vinci and a Chinese 
toy of unknown date. The author's concept consisting in pre- 
sentation in one book all types of propulsors in their histori- 
cal development, should be deemed equally original. 

The book is illustrated with a great quantity of figures 
(303 drawings and fotos) and edited in the form of semi- 
-album.The text of the book is supplemented by : 


— Bibliography (78 items) 

— List of Illustrations 

— Index of Subjects 

— Index of Proper Names 

— Index of Watercraft's Names. 


The presented book was preceded 
by another interesting book of the same author, titled : 
„Two centuries of mechanical propulsion of ships", 
edited (in Polish) by Marpress publishing house in 1997, 
and reviewed in Polish Maritime Research No 2/1997. 


Dowty-Hamilton water-jet propulsive system (1963) 
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Influence of gas turbine controller 
adjustment on ship propulsion system 


behavior in rough sea conditions 


Part 2. The simulation investigations 


Marek Dzida 
Gdansk University of Technology 


ABSTRACT 


The paper presents simulation investigations of influence of gas turbine controller adju- 
stment on ship propulsion system operating in heavy sea conditions, based on the model 
presented in Part I of the paper. The ship propulsion system with two gas turbines dri- 
ving — through mechanical gear — ship controllable pitch propeller, described in Part I, 
was used to analysis. The harmonic disturbances due to sea waves with selected frequ- 
encies were used as input function in computing the system s amplitude characteristics 
at different settings of P and PD controllers. 


Key words : ship propulsion systems, gas turbines, automatic control, sea waving 


INTRODUCTION 


The designing process of the control system of the gas tur- 
bine used as ship's prime mover is governed by the general 
principles of ship power plant designing. However the gas tur- 
bine control system has its specific features. These are : tran- 
sient processes of gas temperature before gas turbine, their in- 
fluence on gas turbine power and efficiency, and on its mecha- 
nical and thermal stresses. In order to limit maximum tempe- 
rature values it is necessary to use some limitations which in- 
fluence transient processes in the entire control system. The 
ship propulsion system with gas turbines driving controllable 
pitch propeller is affected not only by the excitations due to 
settings of propeller shaft's angular speed, and of propeller pitch, 
but also a strong influence on operation of the ship propulsion 
system is exerted by sea waves (sea state). The description of 
the problem, the elaborated simulation model together with 
made assumptions, was presented in Part I of the paper. In this 
part (Part IT) the simulation investigations, carried out by me- 
ans ofthe model, and their results are presented and thorough- 
ly discussed. 


STRUCTURE OF THE INVESTIGATED 
SHIP PROPULSION SYSTEM 
WITH GAS TURBINES 


The ship propulsion system with two gas turbines driving — 
— through toothed transmission gear — ship controllable pitch 
propeller (see Fig.1, Part I) was used to analysis. 

In the propulsion system the propeller shaft's angular speed 
controller and the propeller pitch ratio controller was used. 
For each of the engines the separate control loop was applied. 


In the speed controller two additional limiters were applied : 
that of combustion chamber outlet temperature, which is actu- 
ated when a given threshold temperature is exceeded, and that 
of maximum fuel flow to combustion chamber. The controller 
structure was presented in Part I of the paper. 

The gas turbine in the propulsion system in question is 
a two-shaft turbine with separated power turbine, operating in 
simple open cycle. The parameters of the analyzed ship pro- 
pulsion system are given in Tab.1 [1]. 


Taüb.1. Parameters of the ship propulsion system in question 


EE ONERE шю: 
Rated power of ship propulsion system 462.4 


Rated angular and rotating speed (Орт) | rad/s | 2513 
of a turbine (PT) shaft (npr), | rpm | 24000 
Rated | Rated torque of power turbine — of power turbine 
Rated power of power turbine 231.2 


Propeller Number 1 


|— Res — 

пики — | e [ron] 6000 
ушне |» | - [ems] 
ннан | > [m | as | 
ee [vr 


Rated angular and rotating speed (0р), = 41.89 
of a shaft (np). 6.667 
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For the simulation investigations used was the proportional 
controller (P) or — alternatively — the proportional-derivative 
controller (PD), having the following adjustment coefficients: 


kp - gain (amplification) coefficient of controller 
(for P and PD controllers) 
Tp - differential time constant (for PD controller) 


The harmonic disturbances due to sea waves were applied 
to determine — for assumed disturbance frequencies — the am- 
plitude characteristics at different settings ofthe P and PD con- 
trollers. 


HARMONIC DISTURBANCES 
DUE TO SEA WAVES AFFECTING 
SHIP PROPULSION SYSTEM 


The disturbances from the side of sea waves lead to distur- 
bances of resistance torque of rotary system elements. The 
loads were determined as the quantities of the determinate am- 
plitudes dependent on the disturbance amplitude [2]. 

Such disturbance was introduced to the system as a change 
of the propeller advance speed Vp . The speed Ур is calculated 
in accordance with the relevant formula given in Part I, in which 
the speed change due to sea waves has been accounted for. 

The mean water speed of advance onto propeller, cm, with 
accounted for sea wave action, can be determined from the 
relevant relationship given in Part I. From the continuous sea 
wave spectrum only some its definite frequencies Om were se- 
lected. 

In Tab.2 given are the main wave parameters in different 
sea states, on the basis of which values of cm and о have been 
determined. 

In order to investigate the influence of sea conditions on 
the operation of the ship propulsion system with gas turbines 


Dynamic model 
of GTI gas turbine 


Recording of selected operational 
= parameters of GT1 gas turbine 


Dynamic model 
of GT2 gas turbine 


| Recording of selected operational 
| parameters of GT2 gas turbine 


Tab. 2. Main wave parameters in different sea states [2] 


Wind velocit 


Ё 
Меап 
wave period 


21. 


7 


% | 1.06 | 1.40 | 1.83 | 2.48 | 3.06 | 3.64 
[m/s] 


the amplitude characteristics of the system for a harmonic di- 
sturbance from the side of the propeller (sea wave influence) 
were determined. In the calculations the mean ship’s speed Vo 
was assumed constant, and the characteristics were determi- 
ned for two values of the mean advance speed of water flow 
onto propeller, affected by sea waves : 


Water speed 
on wave surface 


* the speed amplitude c,, which corresponds to 3?B sea state 
(Н = 1.25 m;L-—21.4m)andh, = 0.9 m 
* the speed amplitude cm which corresponds to 8°В sea state 
(H= 11m;L=141 m)andh,=0.9m 


In both cases the characteristics were determined in 
function of the disturbance frequency W,, (wave period Tw). 


INVESTIGATIONS OF CONTROLLER 

ADJUSTMENT INFLUENCE ON SHIP 

PROPULSION SYSTEM OPERATION 
IN HEAVY SEAS 


The simulation model (Fig.7) elaborated by means of the 
MATLAB (and SIMULINK package), presented in Part I and 
[1], was used for the investigations in question. 


namic model 


o 
propeller 


of 
transmission gear 


| 


Dynamic model 


Recording of 
operational 
parameters 
of propeller 


Disturbances due 
to sea waves 


Visualization 
of results 


Fig. 7. Simulation model of ship propulsion system with gas turbines and controllable pitch propeller 
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The amplitude characteristics of the propulsion system were 
determined for two types of gas turbine speed controllers, 
P and PD, controlling the propeller angular speed. The charac- 
teristics were calculated for different controller settings and 
different disturbance amplitudes corresponding to the sea wave 
amplitude at 3°B and 8°B sea state. The calculations were 
performed for the constant value of the propeller pitch ratio : 


H/D = 0.7 and the constant linear speed of ship : Vo = idem. 


The increments of any variable of the propulsion system, X, 
were defined as follows : 


Ax =^=®< 
X 


о 
where : 
X, - initial value of a parameter of propulsion system 
X - instantaneous value of the parameter 


The ship propulsion system's 
amplitude characteristics accounting 
for disturbances due to sea waves 


In the calculations of the system's amplitude characteristics 
in question, the following settings were assumed : 


for the proportional controller (P): kp = 5; 10; 15; (20) 
for the proportional-derivative controller (PD) : 
kp=5 Tp=0.ls;ls;3s. 


The disturbance period T,, was assumed in the range from 
0.1 sto 30 s. The disturbance amplitudes (i.e. the mean propel- 
ler speed of advance cm) were chosen as follows : 


Cm = 0.745 m/s - corresponding to the wave of 


L 7 21.4 m, and H = 1.25 m, (3°B sea state) 


and: 
Cm = 3.450 m/s - corresponding to the wave of 


L= 141 m, and H = 11 m, (8°B sea state). 


In Fig.8 presented are the gas turbine power amplitude cha- 
racteristics for the disturbance due to sea waves corresponding 
to 3?B sea state, at three settings of the P controller ofthe pro- 
pulsion system. 


0.09 


Turbine power amplitude AN; [-] 


1 
Sea wave period T, [s] 


Fig. 8. The gas turbine power amplitude characteristics for 3°B sea state 
conditions, at different settings of the P controller with limiters 


With increasing the mean sea wave period Ту, and the con- 
troller amplification coefficient kp the gas turbine power am- 
plitudes ANpT also increase. In Fig.9 presented are the relative 
amplitude characteristics obtained by relating the power am- 
plitude increment ANprT to the disturbance amplitude cm , for 
the same settings of the P controller and two disturbance am- 
plitudes due to sea waves. From the presented diagrams it ap- 
pears that the courses ofthe relative power pulsations are close 


to each other, hence in further considerations only the frequen- 
cy characteristics for the disturbance amplitude due to sea wa- 
ves corresponding to 8°B sea state (Cm = 3.450 m/s) аге taken 
into account. 


0.05 


o 


m 
iil 


A zaii 
a 0.02 » K 
0.01 


Relative turbine power amplitude 
/с 


0.1 1 . 10 40 
Sea wave period T, [s] 
—O— wave height H = 1.25 m corr. to 3°B sea state 
—— wave height H = 11 m corr. to 8°B sea state 


Fig. 9. The relative amplitude characteristics of gas turbine power, 
АМрт/сњ , in function of the sea wave period Ty, , at the constant 
setting kp — 5 of the P controller with limiters 

In Fig.10 the ship propulsion system amplitude characteri- 
stics accounting for wave-induced disturbances, are presented 
for typical parameters ofthe propulsion system. Fig.10a shows 
the amplitude characteristics for the P controller, and Fig.10b 
— for the PD controller. In both controllers the limiters descri- 
bed in Part I (i.e. the maximum fuel flow limiter and the com- 
bustion chamber outlet temperature limiter) were applied. 

Results for the ship propulsion system 
with P controller (Fig.10a) 

For all quantities characterizing the propulsion system the 
increase of the controller's amplification coefficient kp causes 
the increase of amplitude of the relevant quantity. The ampli- 
tudes of the gas turbine power Npr and the fuel flow mf incre- 
ase along with increasing the sea wave period Ty. The maxi- 
mum amplitudes of the combustion chamber outlet gas tempe- 
rature 04 move towards higher values of the mean wave pe- 
riods along with increasing the controller's amplification coef- 
ficient kp; e.g. for the amplification coefficient kp = 5 the ma- 
ximum amplitude occurs at Ty = 2.5 s, whereas for kp = 15 it 
appears at the wave period Ty = 5 s. 

Results for the ship propulsion system 
with PD controller (Fig.10b) 

Atthe constant value ofthe controller amplification coeffi- 
cient kp the amplitudes of each of the power propulsion sys- 
tem quantities increase along with increasing the differential 
time constant Tp. For all the analyzed system parameters, except 
the fuel flow mp and the propeller speed np, their maximum 
amplitudes occur within the interval between the minimum and 
maximum period of the considered sea wave. At the values of 
the controller differential time constant Tp equal to 1s and 3 s, 
the fuel flow characteristics decrease along with increasing the 
sea wave period Ty. At Tp = 0,1s the transient characteristics 
are almost identical with the transient signals of the P control- 
ler having the same value of the amplification coefficient kp. 
For the propeller speed np the amplitude characteristics differ 
only a little both in function of the sea wave period and of the 
controller's differential time constant settings. 


The influence of controller limiters on the ship 
propulsion system's operation in heavy seas 


Fig.11 shows the transient signals of the selected quantities 
of the considered propulsion system in the case of the use of 
the P and PD controllers both with and without limiters during 
operation of the system in heavy sea conditions (8?B sea state). 
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Fig. 10. The amplitude characteristics of the ship propulsion system affected 
by sea wave disturbances of the amplitude corresponding to 3°B sea state; 
a) for the P controller with limiters b) for the PD controller with limiters 
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Fuel mass flow Am, [-] Propeller torque AM, [-] Propeller speed An, [-] 
N | 
Propeller speed An, [-] 
a A о ON һа 


Temperature АӨ, [-] 


Turbine power AN,, [-] 


a) P controller of k, = 15 
Sea wave period T, = 9.5 s 
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b) PD controller of k, = 5 and T, = 3 s 
Sea wave period T, = 4.5 s 
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Fig. 11. Transient signals of selected quantities of the considered ship propulsion system 
affected by sea wave disturbances of the amplitude corresponding to 8°B sea state 


a) for the P controller; 


b) for the PD controller; 1 — without limiters, 2 — with limiters 
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The results for the system with the P controller are given in 
Fig.11a for the sea wave of the mean period Tw = 9.5 s, and in 
Fig.11b - the results for the system with the PD controller and 
the sea wave of Tw = 4.5 s. The controllers with limiters make 
the correct work of the propulsion system possible, not allo- 
wing for exceedance of gas temperature and stresses in the pro- 
pulsion shaft over the permissible limits both in transient and 
steady states. The work of the system in heavy seas without the 
limiters in the control system is inadmissible. 


The ship propulsion system’s operation 
in heavy seas at different types and settings 
of the speed controller 


The transient signals of the propulsion system operating in 
the heavy sea conditions corresponding to 8°B sea state were 
determined for two mean sea wave periods Ty : 4.5 s and 9.5 s, 
with taking into account two types of the turbine speed con- 
troller: P and PD. The transient signals are presented in Fig.12. 


The transients of the propeller shaft speed np for all inve- 
stigated controller settings show similar form and they do not 
much differ to each other. It results from the fact that it is the 
control signal. 


The transients of the propeller torque M p for the P control- 
ler of kp = 5, and the PD controller of kp = 5 and Тр = 0,1 s do 
not much differ to each other as well. The amplitudes of chan- 
ges are the smallest at both sea wave periods Ty : 4.5 s and 9.5 s 
and all considered controller's settings. The largest amplitudes 
of changes appear for the P controller of kp — 15 and the sea 
wave period Tw = 9.5 s. For the period Тү, = 4.5 s the largest 
amplitudes appear for the PD controller of kp = 5 and Tp = 3 s. 


The amplitudes of fuel flow to gas turbine combustion 
chamber, ту, for the P controller of kp = 5 and the PD control- 
ler ofkp 5 and Tp 7 0,1 s, associated with the influence ofthe 
sea wave of Ty, — 4.5 s, are almost identical. For Ty — 9.5 s at 
the same setting of the P controller they are the smallest. For 
the PD controller the amplitudes increase insignificantly in com- 
parison with the case of the P controller. In both cases the dis- 
turbances due to sea waves do not cause any intervention of 
the limiters to the propulsion system. For the remaining consi- 
dered settings of the controllers the amplitudes are much grea- 
ter, and they are the largest at Tw = 9.5 s for the P controller of 
kp = 15, whereas at Tw = 4.5 s - for the PD controller of kp = 5 
and Tp = 3 s. 


The transients of the combustion chamber outlet tempe- 
rature Occ for the system operating in heavy sea conditions 
are similar to those for the fuel flow my. 


The transients of the gas turbine power Npr at the distur- 
bances due to sea waves are almost the same as the transients 
of the propeller torque My. And, for the P controller of kp — 5 
and the PD controller of kp = 5 and Tp = 0.1s at the wave pe- 
riod Ty, = 4.5 s they are almost identical, and the power oscil- 
lations due to sea waves amont to ANpT = (“Мрт – Npto)/Npto = 
= 25 %; whereas at Ty, = 9.5 s - to АМрт = 32 96 for the P 
controller, and to ANpt = 36 96 for the PD controller. For the 
P controller of kp = 15 and the sea wave period Ty — 9.5 s, the 
power oscillations amount to ANpT = 55 % and are the largest 
out of the results for all considered controller settings. At the 
wave period Ty, — 4.5 s the power changes for the PD control- 
ler of kp = 5 and Tp = 3 s аге the largest (ca 60 96). 
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The performed analysis of the ship gas turbine propulsion 
system operating in heavy seas shows that in the considered 
system : 


e for the P controller kp = 5 should be set; and 

e for the PD controller kp = 5 and Tp = 0.1 s 

e however, for the PD controller at the constant value of its 
amplification coefficient and the increasing differential time 
constant the oscillations of the propulsion system parame- 
ters around their steady values set for operation in heavy 
seas, increase 

e for higher sea state the amplitudes of the considered pro- 
pulsion system parameters increase along with increasing 
the amplification coefficient kp of the P controller 

e forthe PD controller the amplitudes increase with sea state 
increasing, however at a given sea state the maximum am- 
plitudes do not occur at the maximum sea wave period as it 
is the case for the P controller. 


INFLUENCE OF CONTROLLER TYPE 
AND SETTINGS ON OPERATION 
OF THE SHIP PROPULSION SYSTEM 
AT STEP CHANGES OF THE PROPELLER 
PITCH RATIO H/D 


The influence of controller’s type and settings on operation 
of the ship propulsion system in question at step changes of the 
propeller pitch ratio H/D was investigated for the case of calm 
sea. The same propulsion system simulation model as that pre- 
sented in Part I, was used for the investigation. 

Fig.13 shows the transient responses of the ship propulsion 
system on step changes of the propeller pitch ratio H/D within 
the range from 0.7 to 1. The controllers of two types, P and PD 
with limiters and different settings were applied. 

The following control performance criteria were applied 
for the analysis of the influence of the controller’s type and 
settings : 


the minimum control time : 
t, = min 
and the integral performance criteria : 


J, =flett)-e,, dt 
0 


J,= few —e,,) dt 


0 
Fig.14 shows the performance criteria in function of con- 
troller settings determined for the selected parameters of the 
ship propulsion system : np, Npr, Occ. 
From the presented diagrams it results that for the investi- 
gated propulsion system at step changes of the ship propeller 
pitch ratio H/D : 


on the basis of the integral performance criterion J : 
the P controller of kp = 11 should be taken as the optimum one, 
on the basis of the criterion J3 : 
also the P controller of the same kp value; 
however no distinct minimum of this integral appears. 


From the criterion of the minimum control time te of the 
propeller speed np it results that for the P controller a change of 
the controller amplification coefficient kp does not much influ- 
ence the control time. For the PD controller the same criterion 
reveals that the minimum of the control time occurs within the 
interval of the differential time constant Tp from 0.1 to Is. 
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Fig. 12. Transient signals of the considered ship propulsion system affected by sea wave disturbances of the amplitude corresponding 
to 8°B sea state: 1 — P controller of kp = 5; 2 – Р controller of kp = 15; 3З – PD controller of kp —5 and Tp = 0.1 s; 
4 — PD controller of kp = 5 and Tp = 3 s Both controllers fitted with the limiters 
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Fig. 13. Transient responses of the considered ship propulsion system on step changes of the propeller pitch ratio H/D : 
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For such quantities as the gas turbine power Мрт and the FINAL REMARKS 
combustion chamber outlet temperature Өсс, the minimum sta- 
bilization time in response to step changes ofthe propeller pitch 
ratio H/D is achieved: O For the investigated ship propulsion system with gas turbi- 

nes driving - through mechanical gear — controllable pitch 

propeller the optimum controller settings for calm sea con- 
ditions are not optimum ones for the system working in 
heavy sea conditions. 
O In the control loop of the propulsion system the maximum 
a) P controller fuel flow limiter and that saveguarding against increase of 
combustion chamber outlet temperature should be applied. 
0.08 O The optimum setting values for the P and PD controllers in 
regard to sea wave disturbances are lower than the opti- 
mum ones regarding to change of ship propeller pitch. 
It seems reasonable to adjust the controller to the lower 
setting values, i.e. those obligatory for the propulsion sys- 
tem working in heavy sea conditions. In this case the oscil- 
lation amplitudes of its particular quantities are much smal- 
ler than those for the controller adjusted to the optimum 
settings in regard to change of pitch of CP propeller. 


NOMENCLATURE 


The presented analysis demonstrated that : 


in the case of the P controller 
for kp value within the interval from 15 to 20 
in the case of the PD controller 
for its optimum settings : kp = 5 and Tp = Is. 


0.06 


0.04 
Cm - mean speed of water flow onto propeller, caused by sea waves 
e -error 
е - error in steady state 
hp - mean depth of immersion of propeller axis 
kp - controller amplification coefficient 
my - fuel mass flow 
t -time 
Tp - differential time constant 
T, - sea wave period 
V, - linear speed of ship 
Ур - propeller advance speed 
5 10 15 20 Occ- combustion chamber outlet gas temperature 
Gain coefficient k, Mm - frequency of sea - wave - induced disturbance (excitation) 
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propulsion system control time t, in function of the controller settings kp 
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Controllers with limiters; calm sea conditions Laboratory test stand of an aircraft gas turbine 
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Research on influence of some 
ship diesel engine malfunctions 
on its exhaust gas toxicity 


Jacek Krzyzanowski, Kazimierz Witkowski 
Gdynia Maritime University 


ABSTRACT 


This paper is a continuation of the previous articles of the authors, published in Polish Maritime Research 
[1, 2], devoted to pollution of the atmosphere due ship diesel engines in operation. In the paper presented 
are results of the investigations carried out in the Gdynia Maritime Academy laboratory with the use of 
a ship one-cylinder diesel engine combusting heavy fuel oil. Two its possible malfunctions : changes of fuel 
injection pressure and changes of injection advance angle were simulated. Both malfunctions, which might 
occur as a result of incorrect engine regulation or wearing or contamination of engine elements, were 
simulated separately. The presented investigations were focused on finding out relationships between the 
assumed malfunctions and exhaust gas content, especially content of nitric oxides (NO,),very toxic compo- 
unds. The paper is ended by several conclusions arising from analysis of the results which have — apart from 
their cognitive merits — also an utilitarian character as they may be put into practice by ship operators, 
provided that some limitations resulting from the specific conditions of the reported laboratory investiga- 
tions are appropriately taken into account. 


Key words : atmosphere protection, emission of nitric oxides, laboratory test 


INTRODUCTION 


Problems of the sea water protection against pollution from 
ships have been dealt by ship designers and operators for more 
than a half of century. In recent years the marine environmen- 
tal protection was extended over the atmosphere. The revelant 
legal acts aiming at lowering emission of toxic exhaust gas 
components, especially nitric oxides (NOx), are in force (App. 
IV to MARPOL Convention, ISO 8178 standard [3]). 

Technical undertakings to cope with the requirements de- 
termined in the standards go in two directions. The first one 
deals with new engines for which the designers and producers 
search for solutions which would make it possible to obtain 
NOx emission low enough to comply with the standards. The 
other trend is connected with the engines in service. During 
their designing no sufficient attention was paid to exhaust gas 
purity. However today they have to be adjusted to comply with 
the being-in-force requirements for lowering nitric oxides (NOx) 
emission. As it was proved in the publications [1], [2], [4], [5], 
[6], [7] a relatively simple way to limit emission of nitric oxi- 
des is to apply appropriate adjustment operations. 

Despite the adjustment of ship diesel engines is made to 
comply with the exhaust gas purity standards, the emission may 
be changed as a result of changing technical state of the engi- 
nes during service. In this paper presented are results of the 
research on influence of some injection system’s malfunctions 
of a ship diesel engine on its exhaust gas content. The investi- 
gations were performed as an active experiment on the engine 
charged with a heavy fuel oil and operating in the laboratory 
conditions. 
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LABORATORY TESTS 


The test object and stand 


The laboratory tests were carried out on the one-cylinder, 
two-stroke, crosshead engine of longitudinal scavening, which 
was loaded by means of the water brake. The engine was char- 
ged with IF 40 heavy fuel oil. The whole stand together with 
its measurement instrumentation was described in [1]. In com- 
parison to the earlier presented stand a basic change consisted 
in adding an Alfa-Laval viscometer which made it possible to 
continuously measure fuel viscosity and stabilize it at a de- 
manded level. 


Scope of the tests 


The test program was prepared for performing series of 
measurements on the engine with two simulated malfunctions 
of the injection system, as follows : 


> changing the fuel injection pressure - to three selected va- 
lues : 180, 220 and 260 bar 

> changing the fuel injection advance angle - to three selec- 
ted values : -10°, -13° and -16° before the piston’s upper 
dead centre (UDC). 


The measurements were performed within the wide range 
of engine load at the permanent rotational speed of 220 rpm. 
For each of the above selected values six measurements were 
performed at the following engine loads : 25%, 40%, 50%, 
60%, 70% and 80%. The engine load was defined as the per- 
centage ratio of a given engine torque and that rated. 


Description of the tests and their results 5 = 

4 = 

The injection system's malfunctions were introduced- du- ^ 4 z 

ring operation of the engine – by means of a supplementary Е 3 = 

scaled instrumentation of the injector and injection pump. © з tz 

Owing to this it was not necessary to stop the engine before Е 5 e 

each successive test cycle. 5 с 

This way of realization of the measurements guaranteed 5, 1 

runing them in steady conditions. The elimination of influence © 100 


of possible disturbances which could arise from multiple star- 
ting and stopping the engine improved accuracy and reliability 
of the obtained results. 
The measurement results are presented 
in Tab.1 and 2 and in Fig.1 and 2. 


Tab.1. Results of analysis of exhaust gas content 
at three different injector opening pressures A 220 bar pressure == Approx. polynomial for 220 bar 


Fuel | Relative 
injection] engine Exhaust gas content 
pressure} load 
| р | MM, | 
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180 bar pressure Approx. polynomial for 180 bar 


X 260 bar pressure === Approx. polynomial for 260 bar 


Fig.1. NO, content in exhaust gas in function 
of three different injector opening pressures 
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advance angle of -10° Approx. polynomial for advance angle of -10° 


A advance angle of -13° === Approx. polynomial for advance angle of -13° 


Tab.2. Results of analysis of exhaust gas content X advance angle of -16° =_—: Approx. polynomial for advance angle of -16° 
at three different advance angles of fuel injection 


Fig.2. NO, content in exhaust gas in function 
of three different advance angles of fuel injection 


Analysis of the results 


The rated value of the injector opening pressure for the 
tested engine amounts to 220 bar. A change of the pressure 
may be caused by an incorrect initial setting of the injector 
(operator's mistake), or its maladjustment during service 
(a change of technical state of injector's spring). The chosen 
values differed by 20% of that rated, i.e. they amounted to 
180 and 260 bar. 

A drop of the injector opening pressure makes that the fuel 
Is injected to the cylinder a little earlier and the injection lasts 
longer, and spraying the fuel, especially its first portions, is 
worser. The worsening of the combustion process resulting from 
that makes NO, level within the whole range of the applied 
engine loads, lower. 

The rise of the injector opening pressure from 220 bar to 
260 bar results in an increase of NO, emission within the who- 
le engine load range. This can be caused due to improving the 
quality of fuel spraying and combustion process, which leads 
to a higher combustion temperature and — in consequence — to 
a higher NOx content in exhaust gas. 
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An incorrect setting of the injector opening pressure, not 
complying with that recommended by the producer, is always 
considered as a malfunction which should be removed as soon 
as possible. However, from the ecological point of view the 
drop of injector opening pressure makes NOx content in exhaust 
gas lower. Hence it may be a simple way which makes it possi- 
ble to so adjust an existing engine as the NOx emission stan- 
dards to be satisfied. However the accompanying consequen- 
ces should be also remembered, namely : increasing the fuel 
consumption and possible rising the thermal load on engine 
combustion chamber elements. 


The malfunction consisting in a change of the injection 
advance angle may result from an incorrect initial adjustment 
of the engine, or from a random change of that quantity during 
the engine’s operation process. The three following values of 
the advance angle were selected : -16°, -13° (rated) and -10°. 


An increase of the injection advance angle unambiguously 
leads to important increments of NO, content in exhaust gas 
within the entire range of the engine’s load, whereas a delay of 
the injection beginning (a decrease of the injection advance an- 
gle) unambiguously makes NO, content in exhaust gas lower. 


However it should be remembered that both an advance 
and delay of fuel injection starting — in relation to the values 
recommended by the producer and set during statical regula- 
tion of the engine — influences not only the exhaust gas content 
but also other important operational parameters of the engine 
by changing combustion process quality. An increase of the 
advance angle causes an increase of the maximum combustion 
pressure Pmax which may lead to mechanical overloading, and 
its drop — to increasing the fuel oil consumption by dropping 
Pmax and increasing exhaust gas temperature. 


The investigations also revealed (Tab.1 and 2) that the si- 
mulated malfunctions made content of carbon oxide (CO) in 
exhaust gas greater. Both an increase and decrease of the injec- 
tor opening pressure and injection advance angle with respect 
to their rated values led to the above mentioned change. 


CONCLUSIONS 


On the basis of the performed tests 
the following general conclusions may be offered : 


> The engine’s malfunction consisting in an increased injec- 
tor opening pressure makes NO, emission to the atmosphere 
greater. 


> The malfunctions consisting in decreasing the fuel injec- 
tion pressure and in delaying the fuel injection beginning 
favourably reduce NO, emission to the atmosphere. How- 
ever it should be remembered that excessive deviations 
from the rated settings recommended by the engine pro- 
ducer may negatively affect its operational cost by incre- 
asing the fuel oil consumption and its durability (cost of 
engine repairs). 


> Occurrence of the considered malfunctions increases the 
exhaust gas toxicity by increasing content of carbone oxi- 
de (CO). 


2 In practical applications all the conclusions should be con- 


sidered with taking into account the limitations associated 
with the used test object and the realized testing program. 


26 POLISH MARITIME RESEARCH, No 1/2004 


BIBLIOGRAPHY 


1. Krzyzanowski J., Witkowski K. : On possible lowering of the 
toxicity of exhaust gas from ship diesel engines by changing 
their control parameters. Polish Maritime Research, No 1 (19), 
March 1999 

2. Krzyzanowski J., Witkowski K. : On effectiveness of lowering 
the toxicity of exhaust gas from ship diesel engine by 
simultaneous changing its two control parameters. Polish 
Maritime Research, No 1 (23), March 2000 

3. ISO 8178 Standard : Reciprocating Internal Combustion 
Engines — Exhaust Emission Measurement. Parts 1 + 9 

4. Witkowski K., Włodarski J. K., Krzyżanowski J.: Atmosphere 
protection against pollution from sea vessels. Joint Proceedings 
No 14, August 2001 

5. Włodarski J. K., Witkowski K., Krzyżanowski J.: Ecological 
problems of heavy fuel combustion in ship diesel engines. 
Statutory research No 62/DS. Gdynia Maritime Academy, 
Gdynia 1998 and 1999 

6. Vogt R.: The MaK emission concept exemplified by the M32. 
MaK Toplaterne — Diesel Engine Journal for our Business 
Friends, No 68, 1994 

7. Vollenweider J.: Exhaust Emission Control of SULZER Marine 
Diesel Engines. New Sulzer Diesel LTD, Wintherthur, 
Switzerland, 1996 


CONTACT WITH THE AUTHORS 
D.Sc. Jacek Krzyzanowski 
D.Sc. Kazimierz Witkowski 
Faculty of Mechanical Engineering 
Gdynia Maritime University 
81-225 Gdynia, POLAND 
Morska 81-87 
e-mail : jacekk@am.gdynia.pl 


Miscellanea 


TOP KORAB 
awards in the academic 
year 2002/2003 


It is already a tradition that the Society of Polish Naval 
Archtects and Marine Engineers KORAB has granted 
awards for the best M.Sc. theses from among those com- 
mended by the Faculty of Ocean Engineering and Ship 
Technology, Gdansk University of Technology. 


In the academic year 2002/2003 
the awardwinnings were : 


* Mr Lech Wykretowicz, M.Sc. for the project оп: "Pre- 
liminary design of diesel power plant of a 150 000 dwt 
universal bulk carrier of 15 kn speed" 

The project was elaborated under supervision of Mr. 
Jerzy Jamroz, D.Sc., the Department of Ship Power 
Plants. 

* Mr Krzysztof Królak, M.Sc. for the project on : "Pre- 
liminary design of a 2000/1500 TEU cellular conta- 
inership for West-american shipping line" 

The project was elaborated under supervision of Mr. 
Bogustaw Oleksiewicz, D.Sc., the Department of De- 
signing Ships and Ocean Engineering Objects. 


UNDERWATER TECHNOLOGY 


Submarine atmosphere 
regeneration 


Ryszard Ktos 
The Naval University of Gdynia 


ABSTRACT 


In the paper was presented introduction of Polish soda lime (produced by the Chemical 
Company Dwory SA) successfully used for filling the carbon dioxide scrubber of a regene- 
ration system built in KOBBEN class submarines for eliminating carbon dioxide from its 
atmosphere without any additional construction changes. Research on the atmosphere re- 
generation was carried out on а KOBBEN class submarine. The initial CO»? concentration 
was achieved by releasing content of two cylinders filled with CO2/N? gas mixture. Fur- 
ther CO? emission adequate to its maximum production by the ship's crew during bre- 


athing, was simulated by means of a constant emission simulator. The time of the preventive operation amoun- 
ted to about 2 hours. The dynamic absorbing capacity of carbon dioxide was validated experimentally. 


Key words : submarine, live support systems (LSS), soda lime 


INTRODUCTION 
Regeneration procedure 


Air regeneration of disabled/distress submarine (DISSUB) 
is necessary due to tactical or emergency reasons when the 
vessel is not able to emerge on the surface and to ventilate its 
atmosphere. Regeneration of DISSUB's atmosphere consists 
in oxygen enriching and excessive carbon dioxide amount re- 
moving. In most situations carbon dioxide elimination is pro- 
vided by means of soda lime which is a chemical compound 
used for removing carbon dioxide from gas mixtures. Soda lime 
consists of calcium hydroxide Ca(OH)2, sodium hydroxide 
NaOH and/or potassium hydroxide KOH. An approximate 
content of NaOH and/or KOH is 2+5%, of the solid phase. 
On submarines soda lime is placed in a separate absorber con- 
nected to an appropriate ventilation system. A suitable reserve 
of soda lime is also saved on board. 

Procedures of using soda lime on the DISSUBs can be 
different, but it is usually accepted that if carbon dioxide con- 
tent reaches 196,9, the absorption of carbon dioxide has to be 
started. At 296,41 CO» content the DISSUB must be ventilated 
or prepared for abandonment. 

Mean oxygen consumption by one human is assumed to be 
ca 25 dm^hour |, and it is generally accepted that the ratio of 
CO» emission and oxygen consumption, per person, (respira- 
tory quotient) amounts to € = 0.8 [2]. It means that СО» emis- 
sion produced by one person during the respiration process, 
amounting to ca 20 dm ош”, is assumed. These values are 
used to calculate the reserve of soda lime for the DISSUB. 


Soda lime qualities 


Soda lime porosity and bulk density are most important at 
designing the carbon dioxide scrubber but these qualities are 
not available from producers. Therefore, it has been necessary 
to determine soda lime porosity and bulk density. 

Bulk density was measured by weighing in air (1000+ 
+100) cm? of the bed. The soda lime bed mass was (801+20) g 
being equivalent to the soda lime bulk density d = (0.80+0.02) 
Кет” as the buoyancy correction could be neglected due to its 
low value of 1g only, i.e.less than the mass measuring error [1]. 


Tab. 1. Basic qualities of granulated soda lime 


Properties covered by the Polish industrial standard : 
Ne ZN-2001/“DWORY "S.A.-70 


Appearance White or white-grey pellets of 
рр the diameter 2.83.5 mm 
Absorbing capacity with the respect Not less than 26% (by mass) 
to carbon dioxide 
Absorbing capacity with г " 


Negative mesh, square mesh "m 
10+18% (by mass 


Drum strength according " 
to the standard PN-81/c-04533 Moe less Die 220b ymiss) 


Additional properties necessary for fore-designing of the soda lime bed 


Bulk density (0.80+0.02) ke-m? 
0.55+0.02 
Hydraulic diameter (3.470+0.003)-10%m 
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Soda lime bed porosity, defined as the total free inter-mo- 
lecular volume to the total volume occupied by a granula- 
ted material, was determined by water flooding the bed of the 
height | = (75+2) mm and diameter D = (90+2) mm, and the 
volume equal to V = (500-5) cm’, consisting 10236 granules. 
The volume of the water necessary for flooding the bed was 
Vw = (27642) cm” which gave the porosity € = (0.55+0.02), 
typical for granulates. 

Such determination of the mean soda-lime hydraulic dia- 
meter was approximate as the soda-lime granules are additio- 
nally twisted (thus measurement results may differ even by 
10+20%). According to Sauter [1], the mean hydraulic diame- 
ter of soda lime determined from an earlier experiment has 
been equal to (3.470+0.003)-10°° т. 


Scrubber foredesign 


The example calculations of the carbon dioxide scrubber 
are presented in Tab.2 on the assumption that : 


* the air, being in laminar flow, passes through the fully fil- 
led scrubber 

** the resistance (consisted of that imposed by scrubber con- 
struction, the sorbent bed resistance and the scrubber body 
resistance) is relatively small in comparison with the sor- 
bent bed resistance 

* the resistance coefficient for granulates (for one granulate 
fraction) can be defined according to Blake-Kozena's defi- 
nition [1] 

* the particle diameter of the granular bed is determined ac- 
cording to the Sauter's definition 

Ж air viscosity is not dependent on temperature and pressure [3] 

* 1 dm? of soda lime approximately bonds 90 dm? of CO, [3]. 


Tab. 2. Calculation of the sorbent bed dimensions [3] 

Data : 

Flow rate of the breathing gas 

Porosity of the bed 
The granule diameter according to Sauter 
Shape factor 
The breathing gas viscosity (air) 
Flow resistance 


Maximum volume of the absorbed CO; 
Calculations : 


1 - height of the bed 
d, — hydraulic diameter 
Ap- bed resistance 


OBJECT 


The submarine soda-lime scrubber of the internal dimen- 
sions of ca 299 x 278 mm and 276 mm height, was filled with 
(18.73+0.02)kg of the soda lime (Tab.1.) The calculated bulk 
density was d= 0.82 kg: dm" which corresponded to the den- 
sity earlier determined in the laboratory. An overall view and 
technical drawing of the canister is given in Fig.1. 
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Fig.1. Soda-lime scrubber for KOBBEN class submarine: 
a) experimental stand and built in submarine 
b) engineering drawing 


The total flow resistance (the channel flow resistance and 
bed flow resistance) in the carbon dioxide scrubber is measu- 
red with the aid of a differential manometer. For the air flow 
rate (132010) ат? тіп’! the bed resistance was Ap = (100+ 
+2) Pa. The shape factor of soda-lime particles, necessary for 
scrubber foredesign, was determined from the measurement 
results and soda-lime qualities (Tab.1), as shown in (Tab.2). 


METHOD 


Carbon dioxide emission 


Research on the atmosphere regeneration was carried out 
on a submarine of KOBBEN class. Gas cylinders filled with 
(3.50+0.0) kg of CO» and (0.25+0.0) kg of N2 were prepared 
for simulating the initial CO2 concentration in the DISSUB 
atmosphere. 

The initial CO» concentration was realized by releasing 
content of the two cylinders into DISSUB's atmosphere (Fig.2). 
After stabilization of CO» concentration (ca 10 min - Fig.3), 
the stable CO» content amounts to (1.162:0.05)9 6,51 CO» (Fig.3). 
Further СО» emission, adequate to its maximum production 
by the ship's crew during breathing, was simulated by means 
ofa special simulator [1]. Its emission amounted to (10.8+0.05) 
dm? of CO» per min. 


Monitoring system 


The submarine was equipped with a portable atmosphere 
monitoring system. The monitoring system consisted of four 
POLYTRON IR СО» gas analyzers produced by Dragerwerk 


AG Libeck, placed at bow, midship, stern and a hydro-cabin. 
The analyzers were connected through a typical personal com- 
puter and Advantech data acquisition modules into a RS 485 
measuring network (Fig.2). Readouts was recorded every 5 s. 
The PC unit connected to the RS485 measuring network was 
capable of monitoring the submarine's atmosphere by means 
of a special computer program. 
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Fig.2. The system of measuring, CO» releasing 
and dosing it into DISSUB atmosphere 


Regeneration of the DISSUB's atmosphere 


The scrubber of the regeneration system (Fig.2) was filled 
with 13,2 kg of soda lime (Tab.1). Regeneration of the DIS- 
SUB's atmosphere was initiated with start of carbon dioxide 
simulator emission. The regeneration was carried out as long 


Carbon dioxide 
content 


De] 


as one of the analyzers did not show stable CO» content of the 
same value as that before starting a regeneration system. The 
time of the preventive operation amounted to about 2 hours. 

If the dynamic absorbing power of 90 dm? of CO» per 
1 dm? of soda lime, is assumed, the soda lime amount placed 
into the canister (13.2 kg corresponds to 16.1dm' of soda lime) 
should absorb about 1450 dm? of CO». The CO» saturation in 
the soda lime was measured; the measurement results are shown 
in Tab.3. 


Tab. 3. Carbon dioxide saturation in the soda lime samples 
located in different layers of canister 


Moisture Content of 
content absorbed CO; 


Upper part of canister — | — 32 | 68 | 


Sample 
location 


CONCLUSIONS 


* The simple and effective foredesign method of a CO, re- 
moving scrubber was proved experimentally. 

¢ The soda lime complying with the Polish industrial stan- 
dard : ZN-2001/"DWORY'" S.A.-70 was successfully used 
for filling the CO, removing scrubber of the atmosphere 
regeneration system installed in a KOBBEN class sub- 
marine. 

* The dynamic absorbing capacity of carbon dioxide was 
90 dm? of СО, per Idm? of soda lime, due to which the 
effective absorbing capacity of the system in question was 
guaranteed. 
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Fig.3. Rate of carbon dioxide elimination by built in KOBBEN class submarine regeneration system 
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= C onference — — 
SILWOJ 2003 


On 22-24 October 2003 
6? Scientific Technical Symposium on : 


Combustion Engines in Military Applications 


was held at Jurata on Hel Peninsula. It was organized by 
Polish Naval University, Gdynia and Military Engineering 
Academy, Warsaw. 

During the Symposium 57 papers were presented : 14 
of which - at 4 plenary sessions and the remaining - du- 
ring 5 topical sessions on : 


Fuels and oils (6 papers) 

Construction and engineering processes (13 papers) 
Drives and controls (7 papers) 

Operation and Diagnostics (13 papers) 

Technology and construction (4 papers) 


The papers dealt with applications of combustion en- 
gines in army, airforce and navy, and most of them con- 
cerned construction and engineering processes as well as 
operation and diagnostics of the engines. Scientific wor- 
kers of Military Engineering Academy, Warsaw, most con- 
tributed in preparation ofthe symposium materials (22 pa- 
pers), 11 papers were prepared by authors from Polish 
Naval University, and 5 by those of Rzeszów University 
of Technology. The remaining papers were elaborated by 
authors from 12 other universities, institutes and research 


centres. 
8 7 
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SemEco 2003 


Under this acronym 2™ series of seminars devoted to 
ecological problems of operation of combustion engines 
and organized by Naval University, Gdynia, was held in 
the last year. The seminars have been one of the activities 
of Maritime Technology Unit, Section of Transport Me- 
ans, Transport Committee, Polish Academy of Sciences; 
and their main organizer was Cmdr., Prof. L. Piaseczny. 


The following papers 
directly concerned ecological problems : 


Problems of application of biofuels in the light of the 
present development state of piston combustion engi- 
nes — by J. Merkisz (Technical University of Poznan) 

Working and ecological features of the engines sup- 
plied with rape oil, rape-oil methyl esters and their 
mixtures with fuel oil – by Z. Szlachta (Technical Uni- 
versity of Kraków) 

Inventory problems of contaminations emitted by com- 
bustion engines — by Z. Chtopek (Warsaw University 
of Technology) 

Influence of partition of fuel dosage injected to diesel 
engine on its cycle parameters and NO, emission 
by I. Kafar (Polish Naval University) 

Ecological aspects of application of controlled timing 
gears — W. Kozaczewski (Industrial Institute of Moto- 
rization, Warsaw). 


The remaining papers dealt with other 
interesting topics regarding combustion engines : 


* Steam-gas engines applicable to torpedos — their con- 
struction, working parameters and operational featu- 
res — by St. Czarnecki (Polish Naval University) 

* Selected problems of modelling supply processes in 
combustion engines — by M. Sobieszczanski (Acade- 
my of Engineering and Humanistics, Bielsko-Biata) 

* Application of connection graphs and state equations 
to modelling energy systems — by M. Cichy (Gdansk 
University of Technology) 

* Modelling the process of charge inlet into engine's cy- 
linder — by M. Lutowicz (Polish Naval University) 

* Identification of the compression process in ship engi- 

ne cylinder on the basis of indication diagram data 

by S. Polanowski (Polish Naval University). 


тщ 


Diesel engine test stand equipped with МЕХА 9300 Horiba 
exhaust gas analyzer, at Polish Naval University, Gdynia 
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Gdansk University of Technology 


Among 10 faculties of the Gdansk University of Technolo- 
gy (GUT) which its academic community currently consisted 
of 17000 students and 2500 staff, one of the oldest is the Facul- 
ty of Ocean Engineering and Ship Technology (FOEST), be- 
ing the only in Poland which has continuously been educating 
engineers, doctors and professors in the field of marine tech- 
nology since 1945. 

During the 55 years of post-war activities, the FOEST 
educated about 5000 B.Sc. and M.Sc. graduates, as well as 
about 260 persons of Ph.D. and D.Sc. degrees. Eleven of the 
Faculty’s professors were granted the dignity of Doctor Ho- 
noris Causa, and 3 of them — chairs of Polish Academy of 
Sciences. 


ORGANIZATION 
The Faculty is organised of the following departments of : 


Mechanics and Strength of Ship Structures 

Underwater Technology, Hydromechanics and Design of 
Ships 

Ship Manufacturing, Materials and Quality Systems 

Ship Power Plants 

Ship and Ocean Structures Deck Equipment and Systems 
Ship Automation and Turbine Propulsion. 


* t+ +++ 


The main facilities of the Faculty аге: 


* The main building of the Faculty including the administra- 
tion offices, two large lecture theatres and a number of other 
lecture rooms, research rooms including the computer faci- 
lities, computer laboratories, technical laboratories 

Ship Hydromechanics Laboratories 

Ship Automation and Turbines Laboratory 

Ship Deck Equipment and Ship Machinery Laboratory 
Ship Research and Training Centre in Науа 

Laboratory of Ship Hull Structures 

Open Air Laboratory for Ship Hydrodynamics in Наза 
Towing Tank and Cavitation Tunnel 

Ship Machinery Laboratory 

Laboratory for Underwater Vehicles and Composite Struc- 
tures 

Laboratory of Material Science. 
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EDUCATIONAL ACTIVITIES 


The below given tables present the educational specialities 
and lines available in the Faculty. They are divided into two 
groups : 10-semester studies for M.Sc. degree and 7-semester 
studies for B.Sc. degree, both finalized by a diploma thesis. 


Lines of M.Sc. studies 


Ship hydromechanics and design 
Technology and material science 
Structural mechanics and 
construction 

Ship Power Plants 

Automation of Ship Power Plants 
Gas and Steam Turbines 

Ship Deck Equipment 


Shipbuilding 
and Ocean 
Engineering 


Machinery, 
Power Plants 
and Deck 
Equipment 


Lines of B.Sc. studies 


Line of studies 


Steel structures 
Composite structures 
Underwater Technology 


e Design, Construction and Operation 
of Power Plants 
e (Саѕ and Steam Turbines 


Management and Marketing in Maritime Econom 


Training for industrial personnel is one of the main area of 
educational activities of the Faculty. For example the Faculty 
owns the technological licence on *Non-hammering straighte- 
ning of welded steel & aluminium thin-walled structures”. 
Special application courses for engineers and workers enga- 
ged in this area are organised each year. 


OTHER ACTIVITIES 


Besides educational and training activities, the Faculty’s 
scientific personnel carries out multilateral cooperation with 
Polish and foreign universities (more than 70 academic cen- 
tres), takes part in activities of different international organiza- 
tions and associations as their members, participates in scienti- 
fic conferences and symposia. 
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Moreover, during the last 5 years the Faculty organised or 
co-organised about 20 international and domestic scientific con- 
ferences and symposia. 

A usual practice of the FOEST is also hosting 
the young scientists and students from different countries. 


RESEARCH ACTIVITIES 


The Faculty conducts a wide range of research in the area of 
naval architecture, marine engineering and ship technology. 


The main research areas are outlined below : 


> Computer aided ship design methods, non-linear optimisa- 
tion methods in waterborne transport and ship design, com- 
putational geometry and geometric modelling of ship hull 
forms, wave resistance theory, ship hull form analysis and 
optimisation methods, CAD models in design of sailing ves- 
sels and yachts, stability safety analysis of sailing vessels. 

> Stochastic analysis of ship safety in sea environment for stan- 
dardisation of ship safety, criteria for standardisation of fre- 
eboard on the basis of prescribed ship safety, prediction and 
analysis of manoeuvring characteristics of ships, numerical 
simulation of ship manoeuvres, analysis of ports and water- 
ways from the point of view of ship manoeuvring safety, me- 
thods for initial ship design from the point of view of survi- 
vability, analysis of the risk of ship capsizing due to parame- 
tric resonance, experimental analysis of stability of fast plan- 
ning mono-hulls, design and analysis of ship propulsors. 

2 Ship mechanics related to a damaged ship, numerical si- 
mulation of a damaged ship behaviour in natural condi- 
tions, safety of ships in the damage conditions, subdivision 
and damage stability of Ro-Ro vessels, subdivision of tan- 
kers from the point of view of environment protection, sta- 
bilisation of ship motion in waves. 

2 Analysis of ultimate and fatigue strength of ship structural 
elements under multi-axial deterministic and probabilistic 
loading, analysis of ship structure vibrations, application 
of numerical methods including the FEM to ship structu- 
res, determination of sea loads on ships and offshore struc- 
tures, theory and modelling of floating objects as complex 
dynamic systems, determination of hydrodynamic loading 
of liquid-filled tanks, analysis and modelling of segmented 
structures with multi-float support. 

2 The mechanisms of initiation and propagation of fatigue 
and brittle failures in ship materials and structures, deve- 
lopment and conditions for complex destruction mecha- 
nisms of ship structures, methods of analysis of durability 
of ship structures, optimisation of ship manufacturing pro- 
cesses, methods of automated dynamic quality control of 
shipyard manufacturing processes, ecological aspects of ma- 
nufacturing processes in the shipbuilding industry, research 
on ship repair and conversion methods with accounting for 
ecological aspects. 

> Optimisation of reliability of ship propulsion and other sys- 
tems, analysis of damages of ship power plant elements, 
requirements for the power, durability, reliability and safe- 
ty indices of ship power plants and their models, compara- 
tive analysis of classical and novel methods of ship power 
plant design with the use of durability, reliability and safe- 
ty probabilistic models, forming of durability, reliability and 
safety of ship power plants at various design stages, ratio- 
nal procedures of designing the ship power plants for high 
level of durability, reliability and safety. 

2 Ship automation, ship power plant control, dynamic diesel 
turbocharger's systems, steam and gas turbine ship propul- 
sion, dynamics of rotors and gas turbine engines, self-exci- 
ted vibrations of turbomachine rotors. 
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The Faculty is a unique centre in Poland 
that is suitably equipped for carrying out : 


» full-scale tests of hull structure parts under ultimate load 

» research and training in the field of ship manoeuvrability 
by using various large manned models of ships (in coopera- 
tion with the Foundation for Safety of Navigation and Envi- 
ronment Protection in Iława, co-founded by the Faculty) 

> educational projects offerred to the European community 
of shipbuilding students, to be performed in the own rese- 
arch facility located at the Jeziorak Lake in Iława; the sco- 
pe of investigations covers: manoeuvrability, resistance, ship 
stability (conducted for the last 4 years with constantly im- 
proving results) 

» applications of underwater techniques (remotely operated 
vehicles, techniques for naval operations) 

> experiments with the use of an air model turbine rig, and of 
gas turbine test stand. 


INTERNATIONAL COOPERATION 
The Faculty is a member of the following organizations : 


A European Association of Universities in Marine Technology 
A The Association of Polish Maritime Industries 
A Association of Shipbuilding Faculties of Nordic Countries. 


Actually, the Faculty participates in the three EU-funded 
projects, and the Leonardo da Vinci Programme : 


з SANDWICH — Advanced Composite Steel Structures 

з HARDER -Harmonisation of Rules and Design Rationale 

** EURO-MTEC - European Marine Technology Education 
Consortium. 


Tempus, Copernicus, Phare ACE, Tessa and SOCRATES 
ERASMUS are other international or EU programmes 
that the Faculty has already participated in. 


The members of the Faculty staff 
participate in activities of : 


* ISSC — International Ship and Offshore Structures 
Congress 

* ITTC — International Towing Tank Conference 

+ PRADS  - Practical Design of Ship Structures 

* STG — Schiffbautechnische Gesellschaft 

+ GL Poland — Polish Committee of the Germanischer 
Lloyd 

* DNV — International Committee of Det Norske 
Veritas 

* IMO — Committees of the International Maritime 
Organization 

* PAN — Polish Academy of Sciences (7 Committe- 


es and Sections) 


Representatives of the Faculty also act as the members 
of 7 Scientific Societies and Technical Councils 
of different Polish organizations, as well as of : 


European Association of Reliability 

ASME ( USA) 

Association Internationale Cybernetique (Belgium) 
UNESCO - UNISPAR - (University-Industry-Science-Part- 
nership United Nations Educational, Scientific and Cultu- 
ral Organization). 


